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ABSTRACT 

We present a sub-100 pc-scale analysis of the CO molecular gas emission and kinematics 
of the gravitational lens system SDR81 at redshift 3.042 using Atacama Large Millime¬ 
tre/submillimetre Array (ALMA) science verification data and a visibility-plane lens recon¬ 
struction technique. We find clear evidence for an excitation dependent structure in the un- 
lensed molecular gas distribution, with emission in CO (5-4) being significantly more diffuse 
and structured than in CO (8-7). The intrinsic line luminosity ratio is r 8 - 7 / 5-4 = 0.30 + 0.04, 
which is consistent with other low-excitation starbursts at z ~ 3. An analysis of the velocity- 
fields shows evidence for a star-forming disk with multiple velocity components that is consis¬ 
tent with a merger/post-coalescence merger scenario, and a dynamical mass of M(< 1.56 kpc) 
= 1.6 + 0.6 X 10^^ Mq. Source reconstructions from ALMA and the Hubble Space Telescope 
show that the stellar component is offset from the molecular gas and dust components. To¬ 
gether with Karl G. Jansky Very Large Array CO (1-0) data, they provide corroborative evi¬ 
dence for a complex ~2 kpc-scale starburst that is embedded within a larger ~15 kpc structure. 
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1 INTRODUCTION 

Recent large-area surveys have revealed a population of dusty sub¬ 
millimetre (sub-mm) galaxies at high redshift (e.g. Greve et al. 
2004; WeiB et al. 2009) that are believed to be the progenitors 
of present-day massive early-type galaxies (Swinbank et al. 2006). 
These galaxies are characterized by their extreme far-infrared (FIR) 
luminosities, considerable stellar masses and large star-formation 
rates (e.g. Coppin et al. 2008). On the basis of morphological and 
kinematic studies at ultraviolet (UV) and optical wavelengths, it 
has been suggested that most of these galaxies are undergoing a 
major merger event (e.g. Swinbank et al. 2006; Engel et al. 2010). 
In this scenario, their extreme luminosities are the result of merger- 
induced star formation. A competing mechanism, in which star for¬ 
mation is predominantly driven by in-falling gas-rich material has 
also been proposed (e.g. Dekel et al. 2009a; Carilli et al. 2010). 

By using high angular-resolution observations of the molecu¬ 
lar gas component of these galaxies, it is possible to constrain their 
kinematic and morphological properties, and establish the origin of 
their intense star-formation in a way that is neither affected by dust 
extinction nor biased by outflows (Tacconi et al. 2008; Engel et al. 
2010; Karim et al. 2013). While most of these studies have so far 
been hampered by the limited spatial resolution of the observations. 
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the advent of the Atacama Large Millimetre/submillimetre Array 
(ALMA) has made it possible to make sensitive studies of sub-mm 
galaxies on arcsec-scales (Karim et al. 2013). The recent addition 
of the long baseline capacity of ALMA has further improved the 
imaging quality (30 to 300 mas resolution; e.g. ALMA Partnership 
2015a). At the same time, strong gravitational lensing has allowed 
the study of high-redshift sub-mm galaxies on 50-100 pc-scales, 
which is comparable to what is achieved at lower redshift (Swin¬ 
bank et al. 2010). 

In this Letter, we take advantage of the enhanced spatial reso¬ 
lution provided by the long baseline capability of ALMA, in combi¬ 
nation with the magnifying power of strong gravitational lensing to 
study the molecular gas distribution and kinematics of the starburst 
galaxy SDP.81 (H-ATLAS J090311.6+003906; z = 3.042 + 0.001; 
Negrello et al. 2010, 2014; Erayer et al. 2011; Bussmann et al. 
2013; Dye et al. 2014). In a previous paper (Rybak et al. 2015; 
referred to as Paper I), we presented an analysis of the gravitational 
lens mass model and a study of the reconstructed heated dust con¬ 
tinuum emission at a resolution better than 50 pc. Here, we focus 
on a detailed analysis of the CO emission line data that were also 
taken during the science verification of the ALMA long baseline 
array. These data are coupled with high angular resolution imaging 
from the Hubble Space Telescope {HST) and Karl G. Jansky Very 
Large Array (VLA) to investigate the conditions that produced the 
intense burst of star-formation within SDP.81. 
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Figure 1. The image-plane HST (WFC3 F160W) rest-frame UV/optical 
emission from SDR 81 (lensing galaxy subtracted) with the ALMA dust- 
continuum (white contours) and VLA CO (1-0) emission (black contours; 
at 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 relative to the peak emission). For the CO 
(5-4) and (8-7) molecular gas distributions, see ALMA Partnership (2015b). 

2 OBSERVATIONS & LENS MODELLING 

2.1 ALMA CO emission line imaging 

The ALMA spectral line observations of SDR81 and the image pro¬ 
cessing are discussed in detail by the ALMA Partnership (2015b). 
To summarize, the CO (5-4), CO (8-7), CO (10-9) and H 2 O (2o2- 
lii) transitions were observed for between 4.4 and 5.9 h in Bands 
4, 6 and 7. The array consisted of between 22 and 36 antennas with 
baselines ranging from 15 m to 15 km. During our analysis, we 
found that the CO (10-9) and H 2 O emission was too weak to pro¬ 
vide reliable velocity and velocity dispersion maps, so we do not 
consider them further here, and instead we focus on the CO (5-4) 
and CO (8-7) transitions. 

The gravitational lens mass model for SDR81 was derived in 
Paper I by applying a Bayesian pixelated modelling technique (Ry¬ 
bak et al., in prep.) to the visibility data of the Band 6 and 7 contin¬ 
uum emission. This lens model is used for the analysis of the CO 
emission line data here, that is, the lens parameters are kept fixed at 
the most probable values derived in Paper I, while the emission line 
structure and its regularization level are left as free parameters. In 
order to model the two-dimensional spectral maps of the CO lines 
we first subtracted the continuum emission from the wv-dataset us¬ 
ing the line-free channels. As the emission lines are quite broad, we 
averaged the data by 2 channels resulting in a velocity resolution 
lines of 41 km s"^ for both lines. This was found to be a reasonable 
compromise between velocity space sampling (to study the struc¬ 
ture of the lines) and signal-to-noise ratio (SNR) for a robust source 
reconstruction. Also, to limit the number of visibilities used in our 
analysis and to ensure that we compare the various transitions on 
the same angular-scale, a 1 MT cut in the wv-distance was imposed. 
For a Briggs weighting scheme with Robust = 0, this corresponds 
to a synthesised beam with a full width at half-maximum (FWHM) 
of 180 X 131 mas at a position angle of 77.7 deg. 

2.2 Archival 7/57-WFC3 infrared imaging 

SDR81 was observed with the Wide Field Camera 3 aboard the 
HST on 2011 April 4 (GO: 12194; PI: Negrello). The data were 
taken through the FllOW and F160W filters, with exposure times 


of 712 and 4418 s, respectively (see Negrello et al. 2014 for fur¬ 
ther details). The HST imaging detected a single dominant lensing 
galaxy, two gravitational arcs from the background source and ex¬ 
tended low-surface brightness structures around the system. 

The data were reprocessed using the standard astrodrizzle 
pipeline (using a pixel size of 65 mas) and the subtraction of the 
lensing galaxy emission was done using galfit (Peng et al. 2010). 
The image-plane structure of the HST and ALMA continuum emis¬ 
sion (the latter from Paper I) are shown in Fig. 1. Although the 
arcs are offset, due to their different positions in the source-plane 
(see below), the separation between the main and counter-arc is the 
same in both the infrared and mm components, which demonstrates 
that they are at the same redshift. We model the infrared emission 
using the pixelated gravitational lensing code of Vegetti & Koop- 
mans (2009). Consistent with previous studies (see Dye et al. 2014, 
2015; Tamura et al. 2015; Wong et al. 2015), we find that the in¬ 
frared emission has a magnification of yUtot = 11.3 + 0.1. We discuss 
the source-plane properties of the rest-frame UV/optical emission, 
with respect to the dust and gas components, in Section 4. 

2.3 Archival VLA imaging 

Spectral line imaging of the CO (1-0) emission from SDR81 was 
carried out using the VLA in D-configuration on 2010 July 18 
(PI: Ivison; 3.3 X 2.3 arcsec beam-size at position angle -10.1 
degr). The integration time was 1.5 h and 16 VLA antennas that 
had been upgraded with the new receiver systems were used. Fur¬ 
ther details about these observations are describe by Valtchanov et 
al. (2011). We obtained these data from the VLA archive and re¬ 
processed them in the standard manner using the Common Astron¬ 
omy Software Applications (casa) package. In Fig. 1, we show the 
image-plane emission of the CO (1-0) transition, with respect to the 
infrared and mm-continuum emission (see above and Paper I). It is 
clear that the CO (1-0) emission is co-located with the infrared-arcs 
and low-surface brightness infrared emission to the south, demon¬ 
strating that these extended infrared components are at the same 
redshift as SDP.81. The SNR of the VLA dataset is too low to pro¬ 
vide a reliable source reconstruction using our visibility lens mod¬ 
elling technique or to determine if the infrared component to the 
north of SDP.81 is also associated with the lens system. 


3 RECONSTRUCTIONS & ANALYSIS 

The source-plane velocity integrated emission of the CO (5-4) and 
CO (8-7) transitions from ALMA are presented in Fig. 2. These 
maps show clear excitation dependent structure in the molecular 
gas distribution of SDP.81. The CO (5-4) emission has both dif¬ 
fuse and multiple compact structures that extend over ~3 kpc in 
the north-south direction, with an additional low surface brightness 
component that is ~l-2 kpc in size towards the north-east of the 
source. The emission from the CO (8-7) transition is much more 
compact, with an extent of about ~1.5 kpc, and there is only a sin¬ 
gle brightness clump. There is also evidence of a faint tail of emis¬ 
sion that extends by ~2 kpc to the south-west; this tail is not seen 
in the CO (5-4) transition, which implies there is higher excitation 
gas in that region, possibly due to shocks. 

As a check, we compare our reconstructions with the image- 
plane zeroth-moment maps of the counter-arc (ALMA Partnership 
2015b); the counter-arc does not cross a critical curve (see Paper I) 
and so the magnification is roughly constant across this image. We 
find that our source-plane reconstructions are in good agreement 
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Figure 2. The zeroth-moment (left) and uncertainty (right) maps for the 
(upper) CO (5-4) and (lower) CO (8-7) emission lines. The colour-bar is in 
units of Jy kms“^ kpc“^. The lens centre is at co-ordinate (0, 0) and the 
source-plane caustic is shown by the (left) white and (right) grey lines. The 
contours are drawn at 0.2, 0.4, 0.6 and 0.8 relative to the peak emission. 


with the observed image-plane structure of the counter-arc. The to¬ 
tal magnification of the CO (5-4) and CO (8-7) gas components are 
17.0 + 0.4 and 16.9 + 1.1, respectively. 

The total velocity integrated emission from the CO (5-4) 
and CO (8-7) transitions has been determined by directly in¬ 
tegrating the zeroth-moment maps over the total extent of the 
reconstructed source. We find the velocity integrated line in¬ 
tensities are 0.54 + 0.06 and 0.45 + 0.03 Jy km s"^ for the CO 
(5-4) and CO (8-7), respectively. These correspond to line lu¬ 
minosities of ^co( 5 - 4 ) “ 0.91 + 0.10x10^^ K kms“^ pc^ and 
^co( 8 - 7 ) “ + 0.02X 10^^ K kms"^ pc^. This gives an intrinsic 

line luminosity ratio of r 8 - 7 / 5-4 = 0.30 + 0.04, which is consistent 
with other starburst galaxies at redshift 2-3 that typically show low 
excitation (e.g. Scott et al. 2011). 

Ideally, we would like to compare with the line luminosity of 
the CO (1-0) emission. However, simply applying our magnifica¬ 
tion of the CO (5-4) or CO (8-7) transitions to the observed CO 
(1-0) spectrum would be incorrect; the CO (1-0) line is offset from 
the dust continuum emission and so does not have the same mag¬ 
nification (see Fig. 1). Also, the CO (1-0) emission from starburst 
galaxies has been found to be extended by up to ~ 16 kpc (FWHM; 
Ivison et al. 2011), and so the magnification of this gas component 
is likely much smaller than for the CO (5-4) or CO (8-7) emission. 

The spatially dependent structure of the CO emission lines ex¬ 
tends over the caustic curves, which explains the complex velocity 
structure that has been observed in the main arc of SDP.81 (see 
ALMA Partnership 2015b); for example, the image-plane line pro¬ 
file of both transitions are highly asymmetric. In Fig. 3, we present 
the intrinsic line profiles of both the CO (5-4) and CO (8-7) tran¬ 
sitions that are corrected for the differential magnification. We find 
that both transitions actually have more symmetric profiles, con¬ 
sisting of a double-horn profile, typical of rotating disks, with ve¬ 
locity peaks that are separated by ~ 290 kms“^ Again, this result 
is consistent with the image-plane line profiles of the counter arc. 
The CO (5-4) and CO (8-7) emission lines are also offset by about 
-80 km s"^ from the systemic velocity, as defined by the CO (1-0) 
measurement (Frayer et al. 201 1). This could be due to there being 



Figure 3. The intrinsic line profiles of the (left) CO (5-4) and (right) CO (8- 
7) transitions, smoothed with a boxcar of 120 kms“^ width, and assuming 
a systemic velocity corresponding to a redshift of 3.042 and using the radio 
definition of the velocity. 


turbulent in/out flowing higher excitation gas, or that the CO (1-0) 
line also undergoes differential magnification. 

The first- and second-moment maps (velocity and velocity dis¬ 
persion) of the CO (5-4) and CO (8-7) transitions are presented 
in Figs. 4 and 5, respectively. The velocity maps for both emis¬ 
sion lines show a common central solid-body-like component in the 
north-south axis. However, the diffuse and extended component of 
the CO (5-4) and the extended tail of the CO (8-7) show a differ¬ 
ent velocity than expected for solid body rotation. Although the CO 
(5-4) extended feature has a slightly higher velocity dispersion than 
the rest of the emission, the velocity dispersion of the CO (8-7) tail 
feature is significantly different; the large ~200 km s"^ velocity dis¬ 
persion is not due to a single broad velocity component, but is due 
to two velocity components along the same line-of-sight. We find 
that, the velocity dispersion of the CO (8-7) is overall higher than 
the CO (5-4) by about 40 km s"^ 

In order to quantify the asymmetry of the reconstructed ve¬ 
locity (dispersion) structure, we model the moment maps using the 
kinemetry method by Krajnovic et al. (2006). This is based on the 
idea that higher-order moments can be described by a series of con¬ 
centric ellipses of increasing major-axis length; the moment profile 
along each ellipse is then decomposed into a Fourier series with a 
finite number of harmonic terms. From this analysis we find that 
both transitions are characterized by a significant level of rotation 
as well as a significant level of asymmetry at large radii. In the cen¬ 
tral regions of both transitions, the first harmonic term, ki, which 
describes the rotational motion, is found to be substantial and the 
ratio ks/ki , which represents higher-order deviations from simple 
rotation, is consistent with zero, as expected for a disk-like rota¬ 
tion. However, in the outer regions, the ki term abruptly drops and 
at the same time the ks/ki term significantly increases, as expected 
for multiple kinematical components or for post-coalescence per¬ 
turbed disks (Krajnovic et al. 2006; Bellocchi et al. 2012, B12). 

Using kinemetry, Shapiro et al. (2008, S08) have devised a 
criterion to assess whether a galaxy is disk or merger dominated, 
based on the total kinematic asymmetry (Kasym) in the dynamical 
state of the molecular gas. This criterion has been subsequently re¬ 
vised by B12, in order to better distinguish between isolated disks 
and post-coalescence mergers on the basis of large asymmetries 
at large radii. We quantify the total kinematic asymmetry to be 
Kasym,S 08 = 0-65 + 0.10 and Kasym,B 12 = 0.63+ 0.14 for the CO (5-4) 
transition and Kasym,so8 = 0.65 + 0.14 and Kasym,Bi 2 = 0.72 + 0.20 
for the CO (8-7) transition (the uncertainties are a combina¬ 
tion of the errors on the Fourier coefficients and the uncertainty 
due to the choice of centre). These values are consistent with a 
merger/post-coalescence merger interpretation according to both 
criteria (Kasym,so8 > 0-50 and Kasym,Bi 2 > 0.15). Note that although 
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Figure 4. The un-smoothed first-moment (upper-left), second-moment 
(upper-right), first-moment uncertainty (middle-left) and second-moment 
uncertainty (middle-right) maps for the CO (5-4) transition. The lower panel 
shows the modelled (left) first- and (right) second-moment maps from our 
kinemetry analysis. The colour-bar is in units of kms“^, and for the first- 
moment maps is shifted by -80 kms“^ from a systemic velocity corre¬ 
sponding to a redshift of 3.042 (see Fig. 3) using the radio definition of the 
velocity. 



Figure 5. The same as for Fig. 4, but for the CO (8-7) transition. 


the classification criteria have been derived for data with angular 
resolutions that are much lower than the data presented here, we 
find that downgrading our resolution does not affect the derived 
^asym valucs. Finally, the visual criteria by Bellocchi et al. (2013) 
would classify SDR81 as a perturbed disk; while the velocity maps 
show a clear disk-like velocity gradient, the velocity dispersion 
map is relatively smooth without an obvious central peak, indicat¬ 
ing that the disk is offset from the actual centre of mass. 

We estimate the dynamical mass of the disk structure by con¬ 
sidering the maximum radius (1.56 + 0.07 kpc), maximum veloc¬ 
ity (150 + 20 km s"^) and the inclination angle (45 + 8 deg) of the 
outermost kinemetry ellipse that is consistent with an ordered ro¬ 
tation. We find an enclosed dynamical mass of M(< 1.56 kpc) = 
1.6 + 0.6 X 10^^ Mq, which although a factor ~2 smaller than that 
determined by Swinbank et al. (2015), is consistent at the 2cr-level. 
However, as the system is far from being in equilibrium, such dy¬ 
namical mass estimates should be interpreted with care. 

We find significant differences in our visibility-plane recon¬ 
structions (this letter; Rybak et al. 2015) and the image-plane re¬ 
constructions of Dye et al. (2015, see also Swinbank et al. 2015). 
Our reconstructions detect more extended structure within the 
molecular gas distribution that are not seen in the image-plane anal¬ 
ysis. Also, the position and intensity of the compact components 
are seen to vary between the two methods (note the compact struc¬ 
ture varies significantly even within the individual image-plane 
analyses of the 1,1.3 and 2 mm continuum data). These differences 
between the visibility- and image-plane results are likely due to the 
choice of weighting and deconvolution errors in the low SNR im¬ 
ages of the image-plane channel data; the wv-taper and weighting 
that were applied to the image-plane data produces a point spread 
function with significant side-lobe structure that can introduce spu¬ 
rious artefacts in the channel images. For these reasons we believe 
that visibility-plane reconstructions are more robust and less bi¬ 
ased, and the interpretation of intensity and velocity maps from 
image-plane reconstructions requires caution. 


4 DISCUSSION & CONCLUSIONS 

The morphological state of SDR81 is shown in Fig. 6, where it is 
apparent that the reconstructed rest-frame UV/optical stellar com¬ 
ponent {HST F160W) and the combined CO (5-4) and CO (8-7) 
molecular gas component (both from this paper), and the heated 
dust from the rest-frame FIR continuum emission presented in Ra- 
per I form a complex structure. Taking the UV/optical stellar emis¬ 
sion first, we see that this component is dominated by an elon¬ 
gated ~4 kpc structure that is oriented by 45 degrees with respect 
to the north-south axis. There is another brightness peak about 4 
kpc south of the main UV/optical peak, and they are connected by 
a low-surface brightness bridge. Note that the diffuse UV/optical 
structure that is also coincident with the CO (1-0) emission extends 
about ~10 kpc south, as shown in Fig. 6. The CO (5-4) and (8-7) 
molecular gas distribution also looks disturbed, with the region to 
the east being much more diffuse, whereas there is a sharp cut-off 
(shown by the small distance between the intensity contours) in the 
gas towards the west, which forms a ridge at the contact point with 
the stellar component. Interestingly, the extended tail that is seen 
in the CO (8-7) molecular gas distribution neatly fits into the area 
between the main UV/optical component and the bridge feature. 

The ongoing star-formation, as traced by the heated dust emis¬ 
sion, is concentrated to a region of about 2 kpc in size, and is co¬ 
incident with the disk-like rotating structure within the inner part 
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Figure 6. A composite source-plane reconstruction of the rest-frame FIR 
continuum dust (green; see Paper I), the combined CO (5-4) and CO (8-7) 
molecular gas (red) and the rest-frame UV/optical emission (grey). 


of the molecular gas distribution. This component also shows ex¬ 
tension to the east and a sharp cut-off to the western edge. We 
note the detection of a compact molecular gas and star-formation 
component (20-60 M© yr"^ kpc"^) near the brightness peak of the 
UV/optical component, suggesting that star-formation is occurring 
in both structures. However, we find no further evidence of either 
heated dust or excited molecular gas emission across the full 15 kpc 
of the optical component, even though there is an established stel¬ 
lar population and a large reservoir of CO (1-0) available for star- 
formation. This suggests that the current burst of star-formation is 
localised to the northern part of the structure. 

The central structure of SDP.81 is somewhat similar to GN20 
(z = 4.05; Hodge et al. 2012, 2015), where the UV/optical emission 
is offset from the molecular gas component, as traced by lower ex¬ 
citation CO (2-1) and the heated dust emission; these offsets are 
not uncommon in high-redshift galaxies and may be attributed to 
dust extinction of the UV/optical stellar emission or to strong neg¬ 
ative stellar feedback together with gas accretion from gas clumps 
or satellites (Maiolino et al. 2015). Like GN20, there is also evi¬ 
dence for a central star-forming disk at the centre of SDP.81. How¬ 
ever, there are also marked differences, which suggest that the star- 
burst in SDP.81 may be within an interacting/merging system. First, 
our kinemetry analysis of the extended molecular gas distribution 
suggests a post-coalescence perturbed disk. Also, the UV/optical 
emission shows multiple compact and diffuse components over an 
elongated ~15 kpc region. The CO (1-0) data in hand are not sen¬ 
sitive enough or have high enough angular resolution to determine 
the relative velocities of these individual components, but it is clear 
from the overall CO (1-0) emission that they form part of the same 
extended structure. 

The merger scenario has been put forward for the triggering 
mechanism behind the strong star-formation observed within sub- 
mm galaxies. This is based on the observed multiple-component 
morphology and disturbed kinematics of unlensed starbursts and 
from simulations (e.g. Engel et al. 2010; Narayanan et al. 2010). 
Indeed, this conclusion was also recently drawn by Messias 
et al. (2014) who studied the source plane properties of H- 
ALTAS J1429-0028. For SDP.81, the situation is still not clear, 
but further observations of the CO (1-0) emission with the VLA at 
higher angular resolution could measure the relative dynamics of 
the various UV/optical components. This would confirm that they 
form a merging system, with large velocity offsets, or if they are 
part of a clumpy structure with a more chaotic velocity distribu¬ 
tion. 
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